Matter-wave localization in a random potential 
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By numerical and variational solution of the Gross-Pitaevskii equation, we studied the localization 
of a noninteracting and weakly-interacting Bose-Einstein condensate (BEG) in a disordered cold 
atom lattice and a speckle potential. In the case of a single BEG fragment, the variational analysis 
produced good results. For a weakly disordered potential, the localized BEGs are found to have 
an exponential tail as in weak Anderson localization. We also investigated the expansion of a 
noninteracting BEG in these potential. We find that the BEG will be locked in an appropriate 
localized state after an initial expansion and will execute breathing oscillation around a mean shape 
when a BEG at equilibrium in a harmonic trap is suddenly released into a disorder potential. 

PACS numbers: 03.75.Nt,64.60.Cn 



I. INTRODUCTION 

Since the experimental realization of Anderson local- 
ization of a Bose-Einstein condensate (BEC) in a disor- 
der potential [l|, 0], this topic has been the subject of 
intense theoretical and experimental activities. Billy et 
al. [l[ observed the exponential tail of the spatial den- 
sity distribution vifhen a ^^Rb BEC was released into a 
one-dimensional (ID) vifaveguide in the presence of a con- 
trolled disorder created by a laser speckle. Roati et al. 
employed a ID quasi-periodic bichromatic optical lattice 
to observe the Anderson localization of a noninteracting 
^^K BEC. Experimental studies employed bichromatic 
optical lattice [2, Q , shaken optical lattice [^ , cold atom 
lattice Q and speckle potentials [l[ . (The effect of a re- 
pulsive atomic interaction on localization has also been 
studied Q.) Theoretical studies employed bichromatic 
optical lattice [S-Q i '^old atom lattice [13, lll| and ran- 
dom p^ - KL4j potentials among others. 

One can produce quasi-periodic [l5| or random |16l4l8l | 
potentials for Anderson localization by optical means. A 
bichromatic optical lattice is realized by a primary lattice 
perturbed by a weak secondary lattice with incommensu- 
rate wavelength [3. Random speckle potentials are pro- 
duced when light is reflected by a rough surface or trans- 
mitted by a diffusive medium [l9| . In addition to optical 
means, disordered or quasi-disordered potentials could 
also be created by using atomic mixtures [lO| or inhomo- 
geneous magnetic fields [20, [2l[ . Gavis and Castin [l^ 
proposed a controlled way of producing a disordered po- 
tential by using a mixture of two different atomic species. 
In this approach, known as the cold atom lattice, the 
atoms of one of the two species, named "scatterers" , are 
randomly trapped in the sites of an optical lattice. If the 
filling factor is much less than unity, only some of the 
sites will be occupied by one atom and the others will be 
empty. The atoms of the second species, that could be 



weakly affected by the presence of the optical lattice and 
denoted as "test particles" , feel the collisional interaction 
with the randomly-distributed scatterers. The collisional 
interaction will act as a random potential for the test 
particles. 

Here, with variational and numerical solution of the 
Gross-Pitaevskii (GP) equation, we investigate the local- 
ization of cigar-shaped BECs in a ID random speckle po- 
tential and in a disordered cold atom lattice. The random 
speckle potential is modeled by taking identical Gaussian 
spikes distributed randomly [i2|. In this case the spikes 
may overlap. The cold atom lattice is modeled by tak- 
ing a periodic distribution of non-overlapping spikes from 
which a small number of spikes are taken out randomly 
[l(j | . We examine separately two regions of the localized 
BECs. The first region corresponds to the center, where 
the density profile of the localized BEC is quite similar 
to a Gaussian shape. In this region, we also use the vari- 
ational approximation for some analytical understanding 
of the localized state. The second region corresponds 
to the tails of the BEC, where we focus on spatially 
extended wave functions with exponential decay corre- 
sponding to weak Anderson localization [ij, [l^. The 
stability of the localized state is also investigated in both 
cases. To the best of our knowledge this is the first de- 
tailed study of Anderson localization in the presence of 
only a disorder speckle or cold atom lattice potential. 
Most of the previous studies [l^, 123 employed such a 
disorder potential superposed on a harmonic potential. 

In Sec. |TT] we present a brief account of the ID GP 
equation, the disordered potentials, and a variational 
analysis of the GP equation. In Sec. IIIII we present 
numerical results using the split-step Fourier spectral 
method and compare these with a variational analysis. 
In Sec. IIVI we present a brief summary. 



II. ANALYTICAL CONSIDERATION 
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We assume that the trapping potential V{x) is disor- 
dered along the longitudinal x direction with a strong 
harmonic trap in transverse directions. Then it is appro- 



priatc to consider a ID reduction of the three-dimensional 
GP equation by freezing the transverse dynamics of the 
BEG to the respective ground states and integrating over 



the transverse variables [23[. The BEG dynamics of N 



atoms can then be described by the following ID equa- 
tion for wave function u = u{x,t) [23l \24 



.du Id'^u . ,2 ^,, , 



dt 



2dx^ 



u, 



(1) 



with nonlinearity g ~ 2aN/a\^ normalization 
/_ |upda; = 1, and a the atomic scattering length. The 
spatial variable x, time t, and energy are expressed in 
transverse harmonic oscillator units a±_ — y^'h/{muj), 
u!~^ and huj, where m is the mass of an atom and to is 
the angular frequency of the transverse trap. 

The trap V{x) can be modeled by a set of S identical 
spikes randomly distributed along the x axis [12| 



V{x) = Vo^v{x-Xj), 



(2) 



where Vq is the strength of the spike, and v(x — Xj) rep- 
resents the potential of a single spike at position Xj . We 
consider Gaussian spikes |12| 



j{x)^{a^/^) ^cxpf-— j, 



(3) 



with normalization J v{x)dx = 1 and width a. The sta- 
tistical average of the disordered potential ^ is [ijl 



{V{x)) 
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V{x)dx 






(4) 



with D the average spacing between spikes and 2L the 
spatial extension of V{x). The auto-correlation function 
for potential V{x) is defined as [l3] 



C{d) = {V{x)V{x + d)) - {V{x))^ 



(5) 



where (•) represents averagin g as in Eq. Q. An average 
spike height Vs is defined by jla | 



Vs^ 



^JjxiVix) 
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1/2 



(6) 



In presence of str ong disorder a multi-fragmented BEG 
state is expected [13[. Under appropriate conditions, 
however, the shape of localized states may be a single 
fragment, so the variational treatment will be useful [25| . 
The variational approach allows to obtain useful rela- 
tions among the localized state parameters. We consider 
the stationary form (^{x) of the localized state given by 
u{x,t) = exp(— i/ii)(/)(x), with /i the chemical potential. 
The real wave function, (j){x), obeys 



where the prime denotes space derivative. We use the 
off-center variational Gaussian ansatz 




exp 



{X- XqY 

2u;2 



(8) 



with w the width, xq the center and M the normalization 
of the localized BEG. The Lagrangian of the system is 



C 



[m0' - (0') V2 - g0V2 - V{x)^^] dx - A^, 
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(10) 



where pj = xq — Xj. The Euler-Lagrangc equation 
dC/dji — yields the normalization M = 1. We use 
it in the following equations. The remaining equations 
dC/dw = dC/dN ~ yield, respectively, 



1 = — ^ + 



2w^ 



27r (0-2 + ^2) 



E 



1 

4^2 



W\2'K 



E^.' 



2^ - 1 I ^.'(11) 



(12) 



j=i 



and determine the width w and the chemical potential /x. 
The corresponding energy E = J_ [(0')^/2 -|- g(j)^ /2 + 
V (x)4>^]dx is given by 



E = 



9 



4^2 



2w^/2■K 






(13) 



III. NUMERICAL RESULTS 

We perform the numerical integration of GP equa- 
tion ([T]) employing the imaginary- or real-time split-step 
Fourier spectral method with space step 0.04, time step 
0.001 and an initial input Gaussian pulse. The time evo- 
lution is continued till convergence. We checked the ac- 
curacy of the results by varying the space and time steps 
and the total number of space and time steps. Although 
we use a time-dependent approach, the localized states 
are stationary. 



A. A speckle potential 

To generate the speckle potential we use a set of 
random numbers with the standard MATLAB function 
RAND [12] . These numbers are then mapped into the in- 
terval [— i,L] by a linear transformation and denote the 



H(t){x) + 4)"{x)/2 - g4r{x) - V{x)(t>{x) = 0, (7) position {xj) of the spikes. We select S = 300, L = 30 
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FIG. 1: (Color online) (a) Typical disordered speckle poten- 
tial from Eqs. © and dS]) with 5 = 300, L = 30, Vb = 1 
and a = 0.1. (b) The auto-correlation function ^ of the po- 
tential (a). The circles (Num) denote numerical results from 
Eq. (O and the solid line (Fit) is the Gaussian fit (|14p with 
l/fl = 4.2872 and an = 0.128 



and a small width a = 0.1. The spikes are described by 
Eq. ([3]) and the speckle potential by Eq. ^. A typical 
speckle potential V{x) for Vg = 1 is plotted in Fig. [I](a). 
The statistical average of Eq. (|3]) is {V) = 5. Using Eq. 
© the average speckle height is Vs = 4.3519, and the 
spatial auto-correlation function ([5]) is shown in Fig. [1] 
(b) where we also show the numerical fit [12| 



C{d)^VieM-d'/'yR) 



(14) 



with amplitude Vr = 4.2872 and correlation length 
an = 0.128. We study localization with potential V{x) 
of Fig. [T] (a) for different Vq. This potential can pro- 
duce a stationary BEC near a; = 0. We also repeated our 
studies with different speckle potentials generated with 
the same method and parameters as V{x) but by differ- 
ent random processes. Similar conclusions are obtained 
in all cases except that the center of the localized states 
could be different. 

We first investigate the central region of the localized 
BEC, and the numerical and variational results are ex- 
hibited in Figs. [l](a) and (b) where the BEC center xq is 
obtained by a Gaussian fitting of the numerical density 
profile. In Fig. [11(a) we present the density profiles |up 
of the localized BEC for different Vq. Good agreement 
between variational and numerical densities is obtained. 
To understand the effect of Vb, the numerical and varia- 
tional widths (upper pannel) and the chemical potential 
fi and Vs/fJ, (lower pannel) of the localized BEC versus Vb 
are plotted in Fig. [21(b) for g = 0. (The numerical width 
is the root mean square size of the BEC.) Figure [2l (b) 
shows that the width decreases and the chemical poten- 
tial increases with the increase of Vq. However, the rate 
of change is smaller when Vq is large enough (for example, 
Vq > 0.4). It means that the kinetic energy dominates in 
the regime. The quantity Vs/fJ- gives a measure of dis- 
order. However, the values of Vs/n cannot be directly 
compared with those in Ref. [12| where a harmonic po- 
tential is used together with the disorder potential, so 
that a smaller disorder may lead to localization because 
of the additional harmonic potential. 

Next we consider the effects of a small nonlinearity g 
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FIG. 2: (Color online) (a) The numerical (Num) and varia- 
tional (Var) densities |up of the localized BEC versus x for 
different Vb. (b) The numerical and variational widths w, 
chemical potential fi and Vs//i of the BEC versus Vq. 
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FIG. 3: (Color online) The effects of g on the center region of 
the BEC. (a) The width and (b) the energy E and chemical 
potential /i versus g. 



on the localized BEC in the case of a stronger disorder 
(viz. Vq = 1). The width, energy and chemical potential 
of the localized BEC versus g are plotted in Figs. [3l (a) 
and (b). The width, energy and chemical potential in- 
crease with increasing nonlinearity. Increasing nonlinear- 
ity means increasing repulsion and hence increased width. 
Increasing nonlinearity also means a larger positive con- 
tribution to the Hamiltonian and hence a larger energy 
and chemical potential. 

So far, we investigated the central region of the local- 
ized BEC. In the presence of strong disorder, we also 
studied the tail region of the BEC but did not find an 
exponential decay. However, one of the more interesting 
issues of localization is in the presence of a weak disorder 
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FIG. 4: (Color online) (a) The numerical density |up and 
its exponential fit (hatched line) versus x for different Vo and 
p = 0. (b) (the upper panel) The numerical width w and the 
localiztion lengths II (left tail) and Ir (right tail) versus the 
nonlinearity g for Vo = 0.2. (the lower panel) The numerical 
chemical potential /x and energy E versus g for Vo = 0.2. 



when the system is localized due to the disorder nature 
of the potential but not due to the strength of the poten- 
tial. In the presence of a weak disorder, Anderson-like 
localization is expected in the tail region of the BEC. 
The weak disorder of a speckle potential can be acquired 
in two ways. One way is to reduce the strength Vo in the 
potential V{x). Another way is to reduce the disorder of 
the spikes' position. (This second case will be discussed 
in the following subsection). Firstly, we study numeri- 
cally the tail region of the BEC in potential V{x) with 
smaller Vq. Figure|4] (a) shows the effects of Vb on the tail 
region of the BEC for 5 = 0. The hatched lines are ex- 
ponential fitting to the tails by_exp(— |a; — x'\/l). Where 
I is the localization length [l|, [l2] and x' is the center of 
the exponential function. The density profile obtained 
numerically is found to have a clean exponential tail for 
smaller Vq. When Vq is large enough, the exponential 
tail fades away and the tail region approximates the cen- 
tral Gaussian distribution. We also find that the tail is 
modulated by the speckle potential if Vo is small enough. 
In general, the center of the exponential function is dif- 
ferent from the center of the localized state. Also, the 
density profile of the tail is asymmetric around the cen- 
tral position because of the asymmetry of the disordered 
potential. Thus, the localization length for the left tail 
[II] is different from that for right tail {Ir). The stronger 




^u 




■ ' /" ' 


■ 1 1 


(b) 




1 


: 


15 


- 


/'■' 


..■"" ■"••...,^ ^^^,.-^."; 


a 10 




*'*'"*-- 


^, 


...--' 


•■_ 






_l/„=0 


■ 


5 


T 






l'o=»-1 


- 










---V„'0.2 


; 











FIG. 5: (Color online) The expansion of a noninteracting BEC 
in the speckle potential, (a) The time evolution for the density 
profile \uf versus x and t of the BEC for Vo = 0.1. (b) The 
time evolution of the width of the BEC for several Vq. 



the disorder is, the smaller the localization length is. For 
Vo < 0.2 or for Vg < 0.86 corresponding to Vsl [i- < 1-5 
(viz. Fig. [11(b)), we have weak Anderson localization for 
g = 0. These values of Vsl \J- are larger than Vs//i < 0.3 
considered in Ref. [IJl- This is reasonable, as in Ref. 
[1^ the authors consider a speckle potential superposed 
on a harmonic potential. 

We study numerically the effect of the small nonlin- 
earity g on the localization of BEC as well. The results 
are presented in Fig. |4](b), where the numerical width w 
(corresponding to the center region) and the localiztion 
lengths II (left tail) and Ir (right tail) versus nonlinearity 
g are shown in the upper panel. When g is smaller, the 
two localization lengths are similar as shown in the upper 
panel in Fig. |3](b). If g is larger, the asymmetry between 
the two tails becomes larger. The localization lengths in- 
crease sharply with g for g > 0.3. The localization is de- 
stroyed and the BEC becomes a multi-fragmented state 
if g > 0.5 and Vo = 0.2. Compared with Fig. E] (b), the 
lower panel in Fig. 2] (b) shows that the energy is smaller 
when Vo is smaller. The result is accordance with Fig. [2] 
(b). 

We also studied the expansion of a noninteracting BEC 
in the speckle potential. We start from a BEC at equilib- 
rium in the harmonic trap. At time t ~ 0, the harmonic 
trap is suddenly changed into a weak speckle potential 
as in the experiment of Billy et al. [ij. The evolution 
for the density profile of the BEC is shown in Fig. [5] (a) 
where the initial expansion and the consequent breathing 
oscillation of the BEC are illustrated. During the evolu- 
tion, the BEC remains localized and oscillates around 
an equilibrium point. To understand this behavior, the 
time evolution of the width of the BEC is plotted in Fig. 
[5] (b) for several Vq corresponding to a weak disorder. 
During the initial expansion, the BEC expands quickly. 
If Vo = 0, the BEC cannot be localized and expands into 
the whole space (see the solid line) . If the potential has 
weak disorder, however, it will be locked in an appropri- 
ate localized state after a certain amount of expansion. 
After this happens, the system executes breathing oscil- 
lation around a mean shape of the localized state and 
the width remains bounded and does not increase indefi- 
nitely with time. The stronger Vq is, the shorter the time 
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FIG. 7: (Color online) (a) The average spike height of the 
disordered potential Vs of Eq. ([6} versus Si . (b) The chemical 
potential of the stationary localized states versus Vs in case 
of four cold atom lattice. 




I 


/ 


■''-'■'■ ':< 


: 


f 


S^ =0 

...S, .5 

S^ .10 


: 


1 1 . . . . 


. , , . 1 


, 



FIG. 6: (Color online) (a) A typical cold atom lattice (upper 
panel), and the numerical density profile of the localized state 
(lower panel) in exponential scale, (b) The numerical (Num) 
and variational (Var) density profiles |itp of the localized BEG 
with disordered potential (a). 



FIG. 8: (Golor online) The expansion of a noninteracting BEG 
in the weak disordered cold atom lattice, (a) The evolution of 
the density profile of the BEG for Si = 5. (b) The evolution 
of the width of the BEG for several Si . 



for the initial expansion is and smaller is the final width. 



B. A disordered cold atom lattice 

Following the model given by Gavish and Castin \1(^ . 
the disordered cold atom lattice is taken as a set of peri- 
odical spikes with a few random vacancies (holes). To nu- 
merically generate such a potential, we consider S = 100 
periodic spikes within x = [—30, 30] defined by Eq. ([3]) 
with cr = 0.1 and Vq ~ 1. Then we take out a small num- 
ber 5*1 of spikes whose positions are obtained with the 
standard MATLAB function RAND. Figure [6](a) shows 
a typical weakly disordered potential with Si = 5 (upper 
panel) , and the numerical density profile of the localized 
BEG in this potential (lower panel) for g = 0. The mod- 
ulation of the exponential tail of the density profile is 
clearly visible. In Fig. [H] (b) we compare the central part 
of the numerical density with the variational density for 
the potential of Fig. [6] (a). 

In Fig. [7] (a) we plot the average spike height Vs of Eq. 
^ versus 5*1 showing a decrease of Vs with Si. To un- 
derstand the effects of the average speckle height on the 
localized BEG, we numerically calculate the chemical po- 
tential of the localized state versus Vs of Eq. ^ for four 
cases. The disordered cold atom lattices corresponding in 
these cases are created by the same parameters but with 



different random processes. Although the disordered po- 
tentials created by the different random processes possess 
different local characteristics, the global effects are quite 
similar. The value and variation of the chemical potential 
of the system arc smaller if Vs is smaller, corresponding 
to strong disorder. In the case of the weak disorder, for 
example Vs > 1.8, the chemical potentials are very sensi- 
tive to the average speckle height. In this case a large Vs 
corresponds to weak disorder in contrast to the speckle 
potential of Fig. [U where the same denotes strong disor- 
der. 

We also study the dynamics of the localized BEG 
trapped in a weak disordered cold atom lattice. An har- 
monically trapped BEG is suddenly released into a weak 
disordered cold atom lattice shown by Fig. [5] (a). The 
evolution of the density profile and width of the BEG 
are plotted in Figs. |8] (a) and (b), respectively. Because 
of the weak disorder, the BEG can remain as a single 
fragment during the evolution (see Fig. [S](a)). If the po- 
tential is periodical (corresponding to 5*1 =0), the BEG 
can not be localized. For a small Si the width executes 
stable breathing oscillation indicating stable localization. 



IV. SUMMARY 

In this paper, using the numerical and variational so- 
lution of the time-dependent GP equation, we studied 



the localization of a noninteracting and weakly interact- 
ing BEC in a disordered potential. We considered two 
models of the disordered potential corresponding to (i) 
the speckle potential ^ \^ and (ii) the disordered cold 
atom lattice [lO, ll2|. These disordered potentials are 
created by a superposition of narrow spikes randomly 
distributed in space. In the case of a single BEC frag- 
ment, we find that the variational analysis is applicable 
and produces results in good agreement with numerical 
simulation. Our investigation shows that the chemical 
potential (energy) of stationary localized BEC is related 
to the average height of the random potential. The larger 
the average height is, the larger the chemical potential 
(energy) is. For a weakly disordered potential, the lo- 
calized BECs are found to have an exponential tail as 
expected in weak Anderson localization [12[. The expo- 
nential tails are asymmetric and are modulated by the 
disordered potential. The effects of nonlinearity on lo- 
calization are investigated carefully. For a sufficiently 
weak nonlinearity, the increase of the repulsive nonlin- 



earity leads to an increase of the localization length of 
the exponential tail and the increase of the chemical po- 
tential (energy) of stationary localized BEC. For larger 
nonlinearity, the localization will be destroyed and the 
BEC will be in a multi-fragmented state. We also in- 
vestigated the expansion of a noninteracting BEC in the 
disordered potential. We find that the BEC will be locked 
in an appropriate localized state after an initial quick ex- 
pansion and will execute breathing oscillation around a 
mean shape of the localized state when a BEC at equi- 
librium in the harmonic trap is suddenly released into a 
potential with weak disorder. 
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